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ABSTRACT 

We study the stellar population history and chemical evolution of the Milky Way in 

■ a hierarchical ACDM model for structure formation. Using a Monte Carlo method 

based on the semi-analytical extended Press & Schechter formalism, we develop a new 
code GAMETE (GAlaxy MErger Tree & Evolution) to reconstruct the merger tree of 
the Galaxy and follow the evolution of gas and stars along the hierarchical tree. Our 
approach allows us to compare the observational properties of the Milky Way with 
model results, exploring different properties of primordial stars, such as their Initial 
Mass Function (IMF) and the critical metallicity for low-mass star formation, Z^r- 

I In particular, by matching our predictions to the metallicity distribution function of 

^ ■ metal-poor stars in the Galactic halo we find that: (i) a strong supernova feedback 

' is required to reproduce the observed properties of the Milky Way; (ii) stars with 

ff^ , [Fe/H]< —2.5 form in halos accreting Galactic Medium (GM) enriched by earlier 

— ' supernova explosions; (iii) the fiducial model {Zcr = 10"'* .^0, mpopiii = 200Mq) 

provides an overall good fit to the MDF, but cannot account for the two HMP stars 
with [Fe/H]< —5; the latter can be accommodated if < lO^^Z© but such model 

■ overpopulates the "metallicity desert", i.e. the range —5.3 < [Fe/H] < —4 in which no 

stars have been detected; (iv) the current non-detection of metal-free stars robustly 

f~| . constrains either Zcr > or the masses of the first stars mpopiii > O.QM©; (v) 

C^h' the statistical impact of truly second generation stars, i.e. stars forming out of gas 

polluted only by metal-free stars, is negligible in current samples; (vi) independently 

■ of Zcr, 60% of metals in the Galactic Medium are ejected through winds by halos 
^ ' with masses M < 6 x 1O^M0, thus showing that low-mass halos are the dominant 
^ , population contributing to cosmic metal enrichment. We discuss the limitations of our 

study and comparison with previous work. 
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1 INTRODUCTION 

The first stars represent the first source of light, metals and 
dust after the Big Bang. Since they form out of gas of pri- 
mordial composition, they are expected to be metal free. 
For this reason, they have been historically called popu- 
lation III stars (Pop III) to make a distinction from the 
observed population II/I stars, which have typical metallic- 
ities Z > - 10"* and Z ~ 0.02, respectively (Bromm 

& Larson 2004). The impact of first stars on the subsequent 
structure formation history is very strong through difi'erent 
feedback mechanisms: UV photons from the first stars are 
likely responsible for the reionization and heating of the in- 
tergalactic medium (IGM); heavy elements are synthesized 
inside such stars and during the early supernova (SN) explo- 
sions dust grains are formed. Dust and metals are released 
in the interstellar medium (ISM) and eventually enrich the 
IGM through SN winds. The efficiency of all these processes 



depends on the properties of the first stars, which therefore 
control the evolution of the ISM, IGM and of the cosmic 
environment in which galaxies form. 

There are still many puzzling questions about the first 
stars, related to the lack of observational evidences of such 
pristine objects. At present, the primordial Initial Mass 
Function (IMF) represents one of the most discussed top- 
ics in this field. Theoretical studies support the idea of 
a top-heavy IMF, biased towards massive and very mas- 
sive stars (Omukai & Nishi 1998; Abel, Bryan & Norman 
2002; Bromm, Coppi & Larson 2002; Ripamonti et al. 2002; 
Nakamura & Umemura 2002; Schneider et al. 2002; Omukai 
& Palla 2003; Bromm & Loeb 2004; O'Shea & Norman 
2006; Yoshida et al. 2006). Indeed, in the absence of metals, 
the collapsing star forming gas is predicted to fragment in 
~ 10^ Mq clumps, progenitors of the stars which will later 
form in their interiors accreting gas on the central proto- 
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stellar core. Because of the poor understanding of the com- 
plex physics which controls the late phases of accretion and 
protostellar feedback, the final stellar masses are still largely 
uncertain but likely to be in the range SOM© — SOOM© (Tan 
& McKoe 2004). 

In contrast with the above speculations, present-day 
star formation (Pop II/I stars) is known to be character- 
ized by a Salpeter IMF in the mass range O.IMq — lOOM©, 
flattening below masses of O.SSM© (Larson 1998). Since the 
characteristic mass of local Pop II/I stars is ^ IM©, if our 
current understanding of primordial star formation is cor- 
rect, there must have been a transition in the properties 
of star-forming regions through cosmic times. Recent the- 
oretical studies suggest that the initial mctallicity of the 
star-forming gas represents the key element controlling this 
transition (Bromm et al. 2001; Omukai et al. 2001, 2005; 
Schneider et al. 2002, 2003, 2006; Bromm & Loeb 2004). 
Following the evolution of protostellar gas clouds with dif- 
ferent values of the initial metallicity and including dust and 
molecules as cooling agents, Schneider et al. (2002-2006) and 
Omukai et al. (2005) find that when the metallicity is in 
the critical range 10"^ < Zcr/ZQ < 10-" there occurs a 
transition in fragmentation scales from ~ lO'^M© to solar 
or subsolar fragments and argue that this is an indication 
for a transition in characteristic stellar masses. The critical 
metallicity value depends on the fraction of metals locked 
in dust grains, which provide an additional cooling channel 
at high densities enabling fragmentation to solar and subso- 
lar clumps (Schneider et al. 2003, 2006; Omukai et al. 2005; 
Tsuribe & Omukai 2006). Thus, according to this scenario 
the onset of low-mass stax formation in the Universe is trig- 
gered by the presence of metals and dust in the parent clouds 
to levels exceeding Zcr- 

In this scenario, stellar archeology of the most metal- 
poor stars represents a promising way to explore primordial 
star formation. Due to their low metallicitics, metal-poor 
stars arc expected to form out of gas enriched by members 
of the first stellar generation. Therefore, despite of the fact 
that metal-poor stars are observed in the Galactic halo or 
in nearby dwarf galaxies, satellites of the Milky Way (MW), 
these objects might well be living fossils of the first star 
formation episodes in the Universe. One of the most impor- 
tant observational constraints is provided by the distribu- 
tion of stellar metallicities in the halo of the Galaxy, the 
so-called Metallicity Distribution Function (MDF). In their 
early studies, Ryan & Norris (1991) and Carney et al. (1996) 
show that this function peaks at [Fe/H]= —1.6* with a tail 
extending down to [Fe/H] « —4. Although the compilation 
of Ryan & Norris (1991) is relatively old and consists only 
of 300 stars, their MDF is consistent with that obtained by 
more recent and larger samples collected by the HK survey 
of Beers et al. (1992) and by the Hamburg/ESO objective 
prism survey (HES, Wisotzki et al. 2000; Christlieb 2003) 
whose analysis is still underway to quantify completeness 
and selection effects (see Beers & Christlieb 2005 for a thor- 
ough review). The joint HK and HES sample includes 2756 
stars with [Fe/H]< —2. It will be further extended by the 



large number of metal-poor stars contained in the most re- 
cent public data release of the Sloan Digital Sky Survey 
(York et al. 2000), and its planned extension which includes 
the program SEGUE, specifically targeted to collect metal- 
poor stars with [Fe/H]< —2 and to constrain the transition 
from the disk population to the halo. 

The HES has recently revealed the existence of two stars 
with [Fe/H]= -5.3^ (HE0107-5240, Christlieb et al. 2002) 
and [Fe/H]= -5.4 (HE1327-2326, Frebcl et al. 2005). Both 
of these stars show a large overabundance of carbon and ni- 
trogen with respect to iron and other heavy elements, sug- 
gesting a similar origin for their abundance pattern. Essen- 
tially two classes of models have been proposed to interpret 
the observed abundances: (i) the stars were born in a pri- 
mordial environment and had their surfaces polluted at a 
later time (Shigeyama, Tsujimoto & Yoshii 2003; Suda et 
al. 2004); (ii) they had been pre-enriched by the yields of 
one or more SNc (Christlieb et al. 2002, 2004; Umcda & 
Nomoto 2003; Bonifacio, Limongi & Chieffi 2003; Schneider 
et al. 2003; Iwamoto et al. 2005). In the first case, these two 
stars would be low-mass Pop III stars, in conflict with the 
proposed critical metallicity scenario. In the second case, in- 
stead, these stars would be second generation stars enriched 
by the ashes of the first SNe. At present, the available evi- 
dence does not allow to exclude any of the above proposed 
scenarios but shows that the obscrvatioual data on metal 
poor stars are starting to challenge theoretical models. 

This challenge is even more serious as it resonates with 
additional constraints on the primordial IMF placed by 
completely independent observations. In particular, the re- 
cent downward revision of Thomson scattering optical depth 
found in the analysis of the Wilkinson Microwave Amsotropy 
Probe (WMAP) three years data (Page et al. 2006) predicts 
a later reionization epoch, in agreement with a universal 
Salpeter IMF (Choudhury & Ferrara 2006; Gnedin & Fan 
2006). The same conclusion is drawn when theoretical mod- 
els (Salvaterra et al. 2006) are confronted with the observed 
anisotropies in the near-infrared background (Kashlinsky et 
al. 2005). These results tend to favor a relatively unevolving 
Salpeter IMF up to the highest redshifts. 

The motivation of the present study was to explore to 
what extent the shape of the low-metallicity tail of the halo 
MDF can provide insights on the nature of the first stars. 
In particular, we aim at addressing the following questions: 

• What are the implications of the observed low-metallicity 
tail of the halo MDF? 

• Does a sharp MDF cutoff put constraints on the primordial 
IMF? 

• What is the relation between the MDF cutoff and the 

critical mctallicity ? 

• What is the origin of the "metallicity desert", i.e. the 
absence of stars in the range —5.3 < [Fe/H] < —4? 

• What is the statistical impact in the observed samples of 

truly second generation stars, i. e. stars enriched only by the 
ashes of the first stellar explosion, which then can be used 
to constrain the nucleosynthetic products and the masses of 
primordial stars? 



* Hereafter we use the notation [Fe/H]= \og{Npg/NH) — 
log{Npf,/Nf{)Q normalized to the solar values of Anders & 
Grevesse (1989). 



t A recent analysis based on the detection of Fe II lines leads to 
a downward revision of the iron abundance of HE0107-5240 to 
[Fe/H] = -5.7 ± 0.2 (Christlieb, Bessell & Eriksson 2006). 
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Several authors have attempted to interpret the im- 
plications of the observed halo MDF by applying different 
methods. Using a semi-analytical model, Hernandez & Fer- 
rara (2001) deduced that a characteristic Pop III mass in- 
creasing towards high rcdshift is required in order to fit the 
\ow-Z tail of the MDF. Prantzos (2003) also pointed out 
that the stellar metallicity distribution depends sensitively 
on whether instantaneous recycling is adopted or relaxed. 
The possible existence of a low-Z cutoff in the MDF has 
been interpreted by Ocy (2003) in the framework of stochas- 
tic chemical enrichment models, as a result of the various 
metal diflFusion/transport /mixing processes at work in the 
Galactic environment. Along similar lines, Karlsson (2006) 
has noticed that the metallicity desert, if confirmed, could be 
used to extract information about the past Galactic star for- 
mation history. Finally, Tumlinson (2006) presented a chem- 
ical evolution model for the MW which follows the star for- 
mation history in a hierarchical galaxy formation scenario. 
Such model accurately reproduces the halo MDF by Ryan 
& Norris (1991). 

The method proposed in this paper is similar in spirit 
to that of Tumlinson; however, it introduces a number of 
novel features, particularly for what concerns the treatment 
of mechanical feedback and model calibration. Moreover, we 
compare our model results with a new MDF determination 
obtained from joint HK and HES data (including 2756 stars 
with [Fe/H]< -2) by Beers & Christlieb (2006). A detailed 
comparison between the two studios is given in Sec. 6. Using 
the newly developed code GAMETE (GAlaxy MErger Tree 
& Evolution) we follow the gradual build-up of the stellar 
population and metal enrichment of the MW along its past 
hierarchical evolution. The code is based on a Monte Carlo 
algorithm, along the linos of Colo ot al. (2000) and Volon- 
teri, Haardt & Madau (2003), standing on the Extended 
Press Schecther (EPS) theory (Bond et al. 1991; Lacey & 
Cole, 1993). We reconstruct the hierarchical evolution of the 
MW from rcdshift z = 20 and follow, jointly, the chemical 
enrichment and history of its stellar population. The chem- 
ical evolution of both proto-Galactic halos and of the gas 
reservoir into which they are embedded, which we call there- 
after the Galactic Medium (CM), is computed by including 
a physically based description of the mechanical feedback 
produced by the energy deposition of massive stars. In ad- 
dition, this allows us to evaluate the typical masses of ha- 
los which mostly contribute to the GM enrichment and to 
understand how the ejection and accretion of enriched gas 
along the hierarchical tree is reflected in the resulting stel- 
lar MDF. To calibrate the two free parameters of the model 
(star formation and wind efficiency) we compare the inferred 
Galaxy average properties with their observed values. The 
model can then be used to predict the shape of the halo 
MDF for different assumed properties of the first stars, i.e. 
the primordial IMF and critical metallicity. 

The paper is organized as follows: in Sec. 2, wo describe 
the algorithm developed to construct the cosmological^ 



t Wc adopt a ACDM cosmological model with h = 0.73, = 
0.24, f^A = 0.72, n^h'^ = 0.02, n = 0.95 and ag = 0.74, consistent 
with the 3-yr WMAP data (Spergel et al. 2006). We adopt the 
fit to the power spectrum proposed by Bardeen et al. (1986) and 
modified by Sugiyama (1995). 



merger tree of the Galaxy and wc show that it successfully 
reproduces the predictions of the EPS theory. In Sec. 3, we 
discuss the implementation of star formation and mechani- 
cal feedback; Sec. 4 discusses the model calibration. In Sec. 
5 we present the results of the model, which are then com- 
pared with the most recent data in Sec. 6. A brief summary 
and discussion, given in Sec. 7, conclude the paper. 



2 BUILDING THE MERGER TREE 

The most popular model for galaxy formation (ACDM) pre- 
dicts that galaxies form through a scries of merging processes 
of lower mass fragments, in a hierarchical tree picture. Wo 
have used a semi-analytical approach based on the Extended 
Press & Schechter (EPS) theory (Bond et al. 1991; Lacey & 
Cole 1993) to reconstruct the hierarchical merger history 
of our Galaxy. Halo mass distributions drawn from Press 
& Schechter formalism (1974) show an excellent agreement 
with those derived from N-body simulations; its main draw- 
back consists in the lack of spatial information. There are 
several alternative methods to construct merger tree algo- 
rithms from EPS (sec Somorvillc & Kolatt 1999 for a criti- 
cal review); following Cole et al. (2000) and Volonteri et al. 
(2003) (herafter VHM), we have developed a binary Monte 
Carlo code with mass accretion. 

The main building block of the algorithm follows from 
equation (2.15) of Lacey & Cole (1993) re-written in terms 
of the progenitor mass M: 



f{M, Mo) dM : 



1 (dc-dco) 



X exp 



(o-M - o-ifo 



)3/2 



dM 



dM. 



(1) 



This equation gives the fraction of mass in a halo of mass 
Mo at rcdshift zq which, at an earlier time z > zo, belongs 
to less massive progenitors having mass in the range M to 
M + dM (see also Cole et al. 2000). The quantity 5c = Sc{z) 
{ScQ = Sc{zo)) is the critical linear overdensity threshold 
for collapse at rcdshift z (zo); (o-f,^^) is the linear r.m.s. 
density fluctuation smoothed with a top-hat filter of mass M 
(Mo). By construction, the integral of eq. (1) from M = to 
M = Mo must give unity, i.e. all halos were in less massive 
fragments at earlier epochs. Multiplying f{M,Mo) dM by 
Mq/M, the mass fraction can be translated into the number 
of halos per unit mass. 



|^(M, Mo) dM=^ /(M, Mo) dM. 



(2) 



Hereafter we drop the argument (M, Mo) to simplify the no- 
tation. Taking the limit z ^ zo of eq. (2), one can obtain an 
expression for the average number of progenitors in the mass 
range M, M -|- dM into which a halo of mass Mo fragments 
considering a step dz back in time (Cole et al. 2000): 



dN 
dM 



dM = 



dSc 



dcjM 



dM 



dM dz, 



(3) 



where M < Mo and z = zo+dz. This equation represents the 
core of our algorithm. We use it to build a binary merger tree 
that starts from the present day Galactic dark matter halo 
and decomposes it into its progenitors, running backward in 
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time up to redshift z = 20. Given the shape of the matter 
power spectrum in ACDM models, the number of halos di- 
verges as the mass M goes to zero. For this reason, a cutoff 
(resolution) mass Mres , which marks the transition between 
progenitors and mass accretion, is required. At any given 
time-step, halos can either (i) lose part of their mass, cor- 
responding to the cumulative fragmentation into halos with 
M < Mres, or (ii) fragment into two progenitor halos and 
lose mass. At each redshift, the mass below the resolution 
limit accounts for the GM in which the halos arc embedded. 
Note that we refer to mass accretion as the process of mass 
loss in the (backward stepping) merger algorithm. 

Once Mres is fixed, one can compute the moan number 
of progenitors in the mass range Mres < M < Mq/2 that a 
halo of mass Mo at zo is fragmented into during a time-step 
dz 

""^ = 1 dM 

and the accreted mass fraction 

Fa = I -JTJTT- dM. (5) 

J Mo dMMo 

Both these quantities depend on dz (see eq. 3), suitably cho- 
sen to prevent multiple fragmentations {Np < 1). Since the 
number of progenitors decreases with dz, binary algorithms 
generally require a high temporal (redshift) resolution. 

At each time step and for each progenitor mass Mo , we 
generate a random number < R < 1 and we compare it 
with the value Np derived by eq. (4). If Np < R, the halo 
does not fragment at this step. However a new halo with 
mass Mo(l — Fa) is produced, to account for the accreted 
matter. If instead Np > R, fragmentation occurs: a new 
random value in the range Mres < M < Mo/2 is drawn from 
the distribution eq. (3) to produce a progenitor of mass M; 
mass conservation sets the mass of the second progenitor to 
Mo(l — -Fa) — M. This procedure, iterated on each progenitor 
halo and at each redshift, provides the whole hierarchical 
tree of any initially selected galaxy, from the present-day up 
to an initial redshift. 

To give a complete picture of the merger tree algorithm 
we need to describe the quantities involved in the equations 
and specify the adopted free parameters. In particular, in 
eq. (3) the expressions 5c{z) and (^c(z()) represent the values 
of the critical density extrapolated to z = 0. Specifically, 
5c{z) — l.QdiQ/ D{z) where D{z) is the linear growth factor 
(Carrol, Press & Turner 1992), 

5n™(£)_/ 1 209 »m(^) , o4/7. xV^.s 

where fim(2) = nmoil + zf 11-^^0 + {1 + zf^o]~'^ ■ For the 
adopted cosrriological parameters the above equation repre- 
sents a good approximation of the linear theory growth fac- 
tor. On the other hand the linear rms density fluctuation in 
eq. (3), smoothed with a top-hat filter, is given by 

a'^(M) = — / P(k)W^{kR)k'^dk, (7) 
2^ Jo 

where the integral is performed in the k Fourier-space, P{k) 

is the cold dark matter power spectrum and 

W{kR) = 47ri?^[(sin(fcil) - kRco8{kR))/{kR)% 



is the top-hat window function. 

We now need to specify the free parameters of the 
model, Mres and dz. A high value of Mres would be re- 
quired to preserve the binarity of the code (eq. 4) and to 
control the computational cost implied by a very high time 
resolution. However, we need Mres to be small enough to re- 
solve the low mass halos presumably hosting the flrst stars. 
In addition, Mres must be redshift dependent, decreasing 
with increasing redshift, to reproduce the EPS predictions 
(VHM). Indeed, hierarchical models predict less massive ha- 
los at higher redshifts. The halo mass which, at a given red- 
shift z corresponds to a virial equilibrium temperature Tyir 
can be approximated as, 

M(T_.)^10«M,(^) (^) (8) 

(for an exact expression see Barkana & Loeb 2001). Given 
the above redshift dependence, taking Mres oc M{Tvir,z) 
appears physically motivated. For reasons that will be exten- 
sively motivated in the next Section, our analysis focuses on 
halos with mass M > M{Tyir = lO^K, z) = M^iz). Thus, we 
have chosen a resolution mass Mres{z) = M4(2)/10 which 
enables to follow the history of all MW progenitors having 
M > Ma{z) up to a = 20. 

The time-step dz is empirically selected in order to ob- 
tain a good agreement between the EPS predictions and 
numerical results. We started by considering 820 time-steps 
logarithmically spaced in expansion factor between z = Q 
and z = 20. However, since the number of halos in the high 
mass range was found to exceed the EPS predictions, we 
have reduced the time-step by a factor 5 within the redshift 
interval S < z <12. 

Finally, we need to specify the assumed MW dark mat- 
ter halo mass. Following Binney & Merrifield (1998), we have 
assumed Mmw = 10^^ Mq. This values likely represents a 
lower limit to total halo mass, being estimated through ob- 
servations of the most distant galactic satellites currently 
detected. 

Comparison between the code results and EPS predic- 
tions is shown in Fig. 1, where the mean number of Milky 
Way halo progenitors is plotted at different redshifts as a 
function of mass resulting from the average of 200 different 
Monte Carlo realizations of the merger tree. At each redshift, 
the agreement with EPS is satisfactory down to the lowest 
halo masses M > M4(z), i.e. for all the objects of interest 
to our study. Note that, as expected for ACDM models, the 
typical population mass decreases with increasing redshift 
and few collapsed halos exist at = 15, most of the mass 
being below the resolution mass and therefore diffused in 
the GM. 



3 STAR FORMATION AND FEEDBACK 

In hierarchical models for galaxy formation the first star- 
forming halos are predicted to collapse at redshift z ~ 20, 
having masses M ~ 10^ M© and virial temperatures T^ir < 
10* K. The neutral gcis in these minihalos cannot cool via 
atomic hydrogen and relies on the presence of molecular hy- 
drogen, H2, to cool and collapse, ultimately forming stars. 
The role of minihalos in cosmic evolution has been the sub- 
ject of a number of studies, suggesting that their star for- 
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log M [M^] 

Figure 1. Number of progenitors of a MW parent halo of mass Mmw = 10^2^/0 as a function of mass at different redshifts. Histosrams 
represent averages over 200 realizations of the merger tree, with errorbars indicating Poissonian error on the counts in each mass bin. 
Solid lines show the predictions of EPS theory; dashed lines indicate the values of M4{z) at the corresponding redshift. 



mation activity is easily suppressed by radiative and me- 
chanical feedback (see Ciardi & Ferrara 2005 for a thorough 
review). Empirical evidence for this suppression has been 
recently provided by the 3-yr WMAP data. The downward 
revision of the measured optical depth to electron scattering 
of T = 0.09 ±0.03 (Page et al. 2006) is consistent with mod- 
els in which ionizing sources populate halos down to a virial 
temperature of Tvir — 10* K, requiring strong suppression of 
star formation in earlier minihalos (Haiman & Bryan 2006) . 

Consistent with these findings, we have tracked the star 
formation history and chemical enrichment of the Galaxy 
down to progenitor halos with masses Mi{z) and we have 
assumed that at the initial redshift of z ~ 20 the gas in these 
halos is still of primordial composition. This is the result of 
the effect that any given source photodissociates molecu- 
lar hydrogen on scales much larger than those affected by 
its metal enrichment (Scannapieco et al. 2006). In any star- 
forming halo and at each time-step, stars are assumed to 
form in a single burst; the total mass of stars formed is pro- 
portional to the available gas mass with a redshift-dependent 
global efficiency f*{z), Af, = f*Mgas- The initial gas con- 
tent of dark matter halos is equal to the universal ratio, 
Mgas ~ {Qb/^m.)M which decreases after each star forma- 
tion event (see eq. 17). 



3.1 Star formation efficiency 

Star formation in gas clouds occurs in a free-fall time t// — 
{3-K /32Gpy^^ where G is the gravitational constant and p 
the total (dark+baryonic) mass density inside the halo (as- 
sumed to be 200 times denser than the background). Since, 
at each redshift, the Monte Carlo time step At (correspond- 
ing to dz) is At <^ tff{z), it is possible to accurately sample 
time variations of the global star formation efficiency using 
the following approximation; 



At{z) 



(9) 



3.2 Pop III and Pop II/I IMFs 

We have assumed that the stellar IMF depends on the ini- 
tial metallicity of star forming clouds, in agreement with 
the critical metallicity scenario. A star forming halo with an 
initial metallicity Z < Zcr will be referred to as a Pop III 
halo^, and it will host Pop III stars with masses within the 
range of pair instability SNe (SN^^, 14OM0 < m < 26OM0). 
Indeed, these violent explosions are predicted to provide the 
dominant contribution to metal enrichment at the lowest 
metallicities, releasing roughly half of their progenitor mass 
in heavy elements and leaving no remnants. This means that 
all the initial progenitor mass is returned to the surround- 
ing medium and that gas pollution by these objects is very 
strong with respect to the contribution of all stars with 
m < 4OM0 enriching the ambient medium through mass 
loss or Type II SN (SNII) explosions. We have assumed a 
reference value of mpopiii = 2OOAf0 for the mass of Pop III 
stars, but we have also explored the implications of adopt- 
ing the two extreme values of I4OM0 and 26OAf0 in the 
Appendix. 

Conversely, if the initial metallicity exceeds the critical 
value, Z > Zcr, the host halo is referred to as a Pop II/I 
halo and the stars are assumed to form according to a Larson 
IMF; 



$(m) = 4~ oc m exp(— mcut/m), 
am 



(10) 



where e,, physically corresponding to the "local" star for- 
mation efficiency, represents a free parameter of the model. 



with X = —1.35, rricut = O.35Af0 and m in the range 
[0.1 — lOO]Af0 (Larson 1998). This expression represents a 
modification of the Salpeter law which reproduces the ob- 
served present-day stellar population for m > 1 Mq; in the 
low-mass limit, in fact, the IMF behavior is still very uncer- 
tain because of the unknown mass-luminosity relation for 
the faintest stars. The Larson IMF matches the Salpeter 
law for m > 1 Mq while the cutoff at m ~ O.35Af0 can 
explain the absence of brown dwarfs in the observed stellar 
population. 



3 We define Pop III as all the stars with Z < Zcr- In addition. 
Pop III stars are assumed to be massive if Zcr > 0, and dis- 
tributed according to a Larson IMF if Zcr = (see Sec. 6.1) 
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3.3 Instantaneous Recycling Approximation 

Very massive Pop III stars are characterized by a fast evolu- 
tion, reaching the end of their main sequence phase in 3 — 5 
Myr. Conversely, the broad mass range which characterizes 
Pop II/I stars implies a wide range of stellar lifetimes, TsI, 
which vary from a few Myr to several Gyr. In our model, we 
have assumed the Instantaneous Recycling Approximation 
(IRA, Tinsley 1980), according to which stars are divided in 
two classes: those which live forever, if their lifetime is longer 
than the time since their formation Tsi > t(0) — t(zform)', and 
those which die instantaneously, eventually leaving a rem- 
nant, if Tsi < t{0) — t{zform)- The transition mass between 
the two possible evolutions, or turn-off mass mi{z), has been 
computed at any considered redshift. All stars having mass 
m < mi{z) represent stellar fossils which can be observed 
today. The turn-off mass is an increasing function of time 
since [t{0) — t{zform.)] — > when Zform — > 0; in this limit of 
course mi — » lOOM© i.e. all the stars are still alive. Using 
the IRA approximation at each time-step, we can compute 
the number of stellar relics per unit stellar mass formed: 



_ Jo.imq ' 

and the equivalent mass fraction in these stars 

r'"'.*^'' m ^(m)drn 
JO.lMp, ^ ' 



fm* — 



(11) 



(12) 



By definition, there are no stellar fossils of Pop III stars in 
the SN-f-i progenitor mass range. 



3.4 Nucleosynthetic products 

Massive stars can lose mass and heavy elements through 
stellar winds and supernova explosions. Using the IRA ap- 
proximation, we can compute the yield, i.e. the mass fraction 
of metals produced per unit stellax mass formed. 



r^"'!*(® mz(m, Z) <I>(m) dm 



f 

Jo 



lOOMg 



m $(m)dm 



(13) 



as well as the returned fraction, or the stellar mass fraction 
returned to the gas through winds and SN explosions: 

plOOMQ 



R = 



Imi(z f " Wrnim)) $(m) dm 



/oTA^®*('^)"^dm 



(14) 



The quantity mz{m,Z) represents the mass of metals pro- 
duced by a star with initial mass m and mctallicity Z, 
and w{m) is the mass of the stellar remnant. Non-rotating 
Pop III stars in the SN^.-^ domain return all their gas 
and metals to the surrounding medium, i.e. R = 1. We 
have used for these stars the results of Heger & Woosley 
(2002). It is interesting to note that although the total metal 
yield is independent of the progenitor mass and equal to 
Y = 0.45, the iron yield strongly depends on mass, being 
Ype = (2.8 X 10"^^ 0.022, 0.45) for m = (140, 200, 260)Afc?. 
As for Pop II/I stars we have used the grid of models 
by van den Hoek & Groenewegen (1997) for intermediate 
(O.9M0 < m < 8Mq) mass stars and Woosley & Weaver 



(1995) for SNII (SM© < m < 4OM0), linearly interpolating 
among grids of different initial metallicity when necessary. 
We have also followed the evolution of individual elements 
relevant to the present study, Fe and O. 

3.5 Mechanical Feedback 

SN explosions may power a wind which, if sufficiently ener- 
getic, may overcome the gravitational pull of the host halo 
leading to expulsion of gas and metals into the surrounding 
GM. This mechanical feedback has important implications 
for the chemical evolution along the merger tree, as the nu- 
cleosynthetic products of the first stars can be efficiently 
ejected out of the shallower potential wells of the first star 
forming halos and then be accreted by neighboring halos, 
thus being incorporated into the next stellar populations. 

To model mechanical feedback, we compare the kinetic 
energy injected by SN-driven winds 

EsN = e.u.NsN{EsN), (15) 
and the binding energy of the host halo with mass M, 

1 + 



1 GM^ , ,„B3 / Ms 

Eb = = 5.45 X 10^^ erg ( ^ 



2 r„ 



5/3 



10 



h~ 



where Ms = M/W^Mq (Barkana & Loeb 1991). In the first 

equation, is wind mechanical efficiency (i.e. the fraction of 
explosion energy converted into kinetic form); it represents 
the second free parameter (the first is e*, see eq. 9) of the 
model. NsN is the number of SNe in the burst and (Esn) is 
the average explosion energy, which we take to be equal to 
2.7 X lO'^^ erg for SN^-, and to 1.2 x 10^^ erg for SNII. We 
assume that ejection takes place when Esn > Eb, and the 
gas is retained otherwise. We compute the ejected fraction 
of gas and metals as 



{Esn - Eb)/{EsN + Eb). 



(16) 



Hence, the ejection fraction is directly proportional to the 
SN energy provided it is larger than the binding energy. Note 
that, according to this simple prescription, gas and metals 
are ejected with same efficiency. This might not necessarily 

be the case, as shown by several authors (Mac Low & Ferrara 
1999; Fujita et al. 2004); we neglect this complication in this 
work. 

Due to mechanical feedback, the mass of gas and metals 
in a halo can decrease substantially. Following a star forma- 
tion burst and mechanical feedback, the mass of gas left 
in the halo, Mgas, which represents the reservoir for sub- 
sequent star formation events, is related to the initial gas 
mass, Mg"s, and the stellar mass, M*, by 

M,as = [M;L -M^ + i^M.](l - a,j) = 

= M;L(l-/*+-R/*)(l-ae.), (17) 
Similarly, the final mass of metals can be written as, 

Mz = [Mi"(l - /.) + FM.](1 - Qe,) = 

= MjL(^m(l - /») + y/»)(l - aei). (18) 



3.6 Metal mixing 

We have reconstructed the star formation and chemical en- 
richment history of our Galaxy applying iteratively eqs. (17)- 
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Figure 2. Global properties of the MW. Squared points are the model results (200 realizations) as a function of the star formation 
efficiency, e*. Error bars are not shown because negligible in the selected scale. Dashed lines represent the (typical) observed values. Left 
panel: no-feedback model (e„ = 0); Right panel: feedback model with e„ = 0.2. 



(18), together with eqs. (11)-(16), along the hierarchical 
merger tree. It is assumed that during a merger event the 
metal and gas content of two distinct progenitor halos are 
perfectly mixed in the ISM of the new recipient halo. Simi- 
larly, metals and gas ejected into the GM are assumed to be 
instantaneously and homogeneously mixed (we refer to this 
approximation as "perfect mixing") with the gas residing in 
that component. The filling factor Q of the metal bubbles 
inside the volume corresponding to the size of the MW halo 
today, gives an estimate of the validity of the latter assump- 
tion. The perfect mixing approximation is verified when 



Q(^) 



> 1 



where Rh is the bubble radius and 



(A(^)) 



Vmw 



VMwiO){l + z) 



(19) 



(20) 



is the average mean halo separation within the proper MW 
volume Vmw{z), having assumed Vmw{0) ~ 1 Mpc; Nh{z) 
is the total number of halos at redshift z, averaged over 
200 realizations of the merger tree. The value of Ri, can be 
estimated from a Sedov- Taylor blastwave solution: 



Rb{z) 



E{z) 1 1/5 2/5 



■{Pb{z))^ 



(21) 



where Eiz) is the energy released by SN explosions within 
each halo, {pb{z)} is the mean GM density and tni^z) is the 
Hubble time. If we assume that i) star forming halos have 
Mh{z) = Mi{z) and ii) E{z) = Eof*{z)Qt/OmM4{z) where 
Eq = 1.636 X 10"^^ erg/ Mq is the Pop II explosion energy per 
unit stellar mass formed, we find that Q > 1 when z < 11. 
Such limit implies that regions with Z = are no longer 
present beyond that epoch. Additional discussion on this 
issue is given in Sec. 5.2 and Sec. 7. 



4 MODEL CALIBRATION 

Although we have kept our model as simple as possible, it 
includes a number of relatively poorly known (albeit im- 
portant) physical processes that need to be empirically cal- 
ibrated. To this aim we have used the observed properties 
of the Milky Way as a benchmark to fix the best values 
of the two model free parameters, e* and e„. In particular, 
we have compared the results of the simulations at redshift 
z = with the following observations: 

• Gas metallicity. The Galactic disk has a mean metal- 
licity Zgas — Zq. As this component contains the majority 
of the MW gas, we use this value as representative of the 
whole system. 

• Stellar metallicity. The typical mean value quoted for 
this quantity is Z* ~ Zq. This value has been derived by 
weighting the mean metallicity of each galactic stellar com- 
ponent with its corresponding mass. 

• Stellar mass. Contributions to the stellar mass come 



from the disk {M^ 



[m: 



rbulge 



(4 - 6) X IO^^Mq), the bulge 



(0.4 - 1) 10'"Mq) and the halo (M," 



(0.2 ~ 1) 10"' Mq) components (Dehennen & Binney, 1998; 
Brown et al. 2005), summing up in a total mass of ~ 
6 X 10^" M0. 

• Gas-to-stellar mass ratio, Mgas/M* = 0.13. The mass 
of gas has been derived using the observed mass of HI and 
HII regions of the Galaxy, Mgas = Mhi + Mrii ~ (6 -|- 2) x 
10^ Mq ~ 8 X 10'' Mq (Stahler & Palla 2004). 

Moreover, if mechanical feedback is considered, two addi- 
tional quantities become meaningful: 

• Baryon to dark matter ratio, /(, = (M* -|- 
Mgas)/MMW = 0.07. 

• GM metallicity. Although the exact value for this quan- 
tity is somewhat uncertain, high-velocity clouds (supposedly 
constituted by leftover gas from the Galactic collapse, and 
currently accreting onto the disk) may be used as a rea- 
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model 




< Z*/Z© > 


ZgmIZq 


Mgas/M^ 




fb 


no feedback: e, = 0.5 , = 


1.261 


1.13 




0.115 


1.71 




feedback: e» = 0.7 , ew = 0.2 


0.93 


1.01 


0.41 


0.11 


0.85 


0.085 


observed values 


1 


1 


0.25 


0.1 


1 


0.066 



Table 1. Averaged values of the global properties of the MW derived using the two best- fit models; no-feedback e« = 0.5, e-w = 0, and 
feedback e* = 0.7, = 0.2. The values have been obtained averaging over 200 realizations of the merger tree. For comparison, we also 
quote the indicative observed values (see text). 



sonable template. For these objects, Ganguly et al. (2005) 
derive [0/H]~ —0.66 ~ 0.25Zq, a value that we take as an 

educated guess for Zgm- 

Obviously, in study cases in which mechanical feedback is 
turned off for comparison sakes, it follows that Zqm = 
and ft — flb/i^m ~ 0.16. In principle, an acceptable fit to 
the above MW properties can be achieved both for models 
including mechanical feedback and neglecting it (see next 
Sections). However, the best values of the free parameters 
are different in the two cases. We will show later how it 
is possible to assess the role of feedback by using the MDF 
properties. Note that the calibration procedure is completely 
independent of Zcr and Pop III IMF. 



4.1 No Feedback Model 

In Fig. 2 we compare the observed global properties of the 
Galaxy with the results of our model as a function of the 
star formation efficiency, e*. We start by analyzing the no- 
foodback case, i.e. when e„ = 0. A good overall agree- 
ment with the four observed quantities is obtained for rel- 
atively high values of e* > 0.1. In fact, the derived values 
of Zgas/ZQ, Z^/Zq, M^/M^^ rapidly increase when the 
star formation efficiency increases from e* = 0.01 to 0.1, 
and remain nearly constant thereafter (note that the value 
of Mgas/Mt. actually decreases with e« in this range). These 
trends can be understood very simply. The gas mass de- 
creases due to star formation (eq. 17), forcing the gas frac- 
tion to decrease with increasing e». Thus, Af, saturates be- 
yond e« > 0.1 because of the limited gas mass available. 
Given the proportionality between Mz and M, (eq. 18), 
(Zt) shows the same behavior. Conversely, the mass of met- 
als in the gas is determined by the balance between met- 
als synthesized and ejected by stars, and metals locked into 
newly formed stars (eq. 18); when the second effect domi- 
nates, it produce a slight decrease of Zgas- The best-fit model 

requires a high star formation efficiency, e* — 0.5^. The de- 
rived values for all the observed properties of the MW are 
listed in Table 1. The values shown represent an average 
over 200 realizations of the merger tree. 



^ Because of the assumed star formation law (eq. 9), e* = 0.5 
is only the instantaneous star formation efficiency and does not 
imply that half of the gas mass is converted into star in a single 
burst. Rather, for this model /«(2) varies in the range (0.06 — 4) X 
10~^ in the considered redshift range. 



4.2 Feedback Model 

In the right panel of Fig. 2, we show the global properties of 
the MW derived using a feedback model with ~ 0.2. For 
low values of e* , the effect of mechanical feedback is negligi- 
ble; similarly to the no feedback model ft, il,b/il,m ~ 0.16, 
and Zgm ~ 0. These results are coherent with our feedback 
prescriptions (eqs. 15 and 16) since the efficiency of mechan- 
ical feedback depends on the number of SN explosions and 
therefore on the level of star formation activity. As expected, 
fb decreases with e» while Zgm increases. Note that in the 
range 0.1 < e, < 1 the feedback model closely resembles the 
no-feedback one, the major difference being the lower sat- 
uration values for the stellar mass and mctallicity. Instead, 
the metallicity of the gas is nearly unaffected because met- 
als and gas are ejected into the GM with the same efficiency 
(see eqs. 17 and 18). The best-fit model requires e, — 0.7 
(we stress again that this is not the global star formation 
efficiency, expressed by f*{z), see eq. 9), and reproduces the 
MW global properties better that the no-feedback model 
(see Table 1). 



5 MODEL RESULTS 

5.1 Metallicity Distribution Function 

After having ensured that our models reproduce the global 
properties of the MW, regardless of the value of Zcr and of 
the PopIII IMF, we are now ready to explore the predictions 
of the calibrated models. The result shown below are the 
result of averages over 200 realizations of the MW merger 
tree; in some cases we use a single representative realization 
for illustration purposes. 

In Fig. 3 we compare the MDF (i.e. the number of MW 
relic stars as a function of their iron abundance, [Fe/H], 
a quantity commonly used as a metallicity tracer) derived 
using the two best-fit models, with and without mechanical 
feedback (see Table 1), for the specific value Zcr = W~*Zq 
and m PopIII = 200Mq. 

The first striking difference is that no stars with 
[Fc/H] ;S — 1.5 arc produced in the no-feedback model. The 
reason is that the gas metallicity is not efSciently diluted in 
closed box models. A simple estimate elucidates this point. 
Following a burst of Pop III stars, the gas metallicity of the 
host halo is raised up to 



Z = 



Mz 



YM, 

Maas 



= r/.(z) > 0.146Ze 



(22) 



where we have used Y = 0.45 and the smallest f*{z) 
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Figure 3. Left panel: MDFs derived using the no- feedback model (shaded histogram), and the fiducial feedback model (unshaded) 
both summarized in Table 1. Histograms have been obtained as average over 200 realizations. Both models assume Zcr = 10~^Zq and 
f^popiii = 200Mq. Right panel: Age-metallicity relation for the halos derived using the fiducial feedback model. The starred (solid) 
points represent Pop III (Pop II/I) star-forming halos. For clarity results are shown for a single realization. 



6.5 X 10~ value corresponding to e, = 0.5 (appropriate 
for the no- feedback case). Note that the final metallicity 
depends only on the metal yield and on the star formation 
efficiency. From the previous equation it is easy to see that 

the resulting abundance is I 



[Fe/H] ~ log {Z/Zq) - 0.257 , 



1.09 



This value is larger than the low [Fe/H] cutoff in the dis- 
tribution shown in Fig. 3. What is the reason for such 
discrepancy? In the absence of mechanical feedback, met- 
als can only be diluted through mass accretion of metal- 
free gas and through mergers with unpolluted progenitors 
(which at most can halve the gas metallicity Z if they have 
a comparable mass). However, dilution purely provided by 
accretion of unpolluted gas is not sufficient to account for 
the low-Z tail of the MDF, which extends to values [Fe/H] 
~ —4, with two outliers at [Fe/H] < —5. The same con- 
clusion is achieved using different mass values for SN-y.^ or 
assuming a Larson IMF for metal free stars though, in both 
cases, the different value of the yields cause a shift in the 
MDF cutoff. For example for the Larson IMF {Z^r = 
model) Y ~ 0.01 and Ype ~ 5 x 10~* and the cutoff moves 
to [Fe/H]~ —2.44 » —4. A similar result is obtained for 
mpopiii — 140 M© shown in the Appendix. 
In conclusion mechanical feedback is required to effectively 
dilute metals inside halos independently of the Zcr value and 
of the assumed PopIII IMF. 

Heavy elements ejected by SN winds mix with the ex- 
ternal GM and can be accreted by other halos, raising their 



II This relation follows from the definition of [Fc/H] that leads 
to [Fe/H]= log(Z/Zo) + l-044-|-log(yFe/yz), where for Pop III 
stars with m = 2OOM0 the iron yield is approximately Yp^ = 

o.osr. 



metallicity. Contrary to self-enrichment by previous star- 
burst within the same host halo and to hierarchical enrich- 
ment by contaminated progenitors, this accretion-driven en- 
richment mode is present only when mechanical feedback 
operates. Accretion of metal-enriched GM onto a halo of 
primordial composition sets the initial conditions for the en- 
vironment in which most metal-poor stars will later form. 
Indeed, the unshaded histogram in Fig. 3 shows that when 
mechanical feedback is active the MDF extends down to 
[Fe/H] ^ log(^cr/^©) as low-mass stellar relics can be pro- 
duced only in gas clouds with Z > Zcr- 

The right panel of Fig. 3 shows how different enrichment 
modes operate along the merger tree in a single representa- 
tive realization. Each point denotes the metallicity of a halo 
at a given redshift; in this plane halos move from left to right 
i. e. they increase their metallicity. The evolution of each ha- 
los is followed along its track until it suffers a merger event. 
Two families of curves can be identified, corresponding to 
self-enrichment and accretion mode, respectively, as marked 
by the labels in the Figure. In general, at the onset of the 
first burst of star formation a halo could be either metal-free 
([Fe/H] = — oo) or it could have a non-zero metallicity if it 
has accreted gas from the GM previously polluted by met- 
als ejected from other halos**. In the first case, the halo will 
host Pop HI stars; if metals produced are retained after their 
explosion the halo becomes self-enriched and moves directly 
on the corresponding self-enriched track. In the second case, 
instead, depending on the relative ratio of the GM metal- 
licity Z and Zcr ~ 1Q~'^Zq the stars form could be either 
Pop HI or Pop II. In this particular realization at z ~ 11 
the GM metallicity crosses the Zcr value and sets the termi- 



** The same conclusion applies to halos that are just about to 
gravitationally collapse and virialize. 
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nation of the Pop III formation epoch in accretion-enriched 
halos. Such halos lie on the curve marked as "accretion". 
However the second star-formation episode will push again 
these halos onto the corresponding self-enriched tracks. The 
accretion-mode is crucial in order to form long-lived stars 
with [Fe/H]< —2.5. The accretion-enriched branch ends at 
2 ~ 1, when merger events dominate the evolution of MW 
progenitors. The last point of the branch corresponds to 
[Fe/H] ~ —1, reflecting the metallicity of the GM at the 
same redshift (see Fig. 4). 

The analysis of this Figure helps interpreting the MDF. 
Our model predicts that stars with [Fe/H]< —2.5 form only 
in halos enriched through the accretion-mode. Conversely, 
iron-rich stars, with [Fe/H]> —1, form in self-enriched halos. 
In the intermediate range, —2.5 <[Fe/H]< —1, stellar relics 
originate from both self-enriched and accretion-enriched ha- 
los, and their contribution to the MDF cannot be separated. 
However, at a given [Fe/H], stellar relics along the accretion- 
enriched branch are systematically younger. 

Note that another possibility to increase metal dilution 
along the merger tree is to relax the hypothesis that all 
Pop III stars end up as S'N^-y {f^-, = 1). Indeed if a fraction 
of Pop III stars forms with m > 260Mq or m < 140Mq 
{fjj < 1) and ends up as black holes (Heger & Woosley 
2002), metal enrichment can be substantially reduced. How- 
ever, Schneider et al. (2006) have recently studied the de- 
pendence of chemical feedback from the parameter f^-, (thus 
from the assumed Pop HI IMF), showing that f^-, = 1 is cur- 
rently favored by observations of the cosmic star formation 
history, high redshift number counts and cosmic reioniza- 
tion. 

In conclusion mechanical feedback is necessary to repro- 
duce both the global properties of the MW and the observed 
range of [Fe/H] in halo stars. Hereafter we will refer to the 
best-fit feedback model of Table 1 as our fiducial model. 



5.2 Metallicity evolution of the Galactic Medium 

In Fig. 4 we show the evolution of GM iron and oxygen abun- 
dances together with the specific contribution to [0/H] by 
Pop HI and Pop II stars. For reference we also plot the ob- 
served [O/H] value measured in high- velocity clouds (Gan- 
guly et al. 2005) which are taken as an indicator of left- 
over GM in the MW halo. These results refer to our fiducial 
model with Zcr = 10~'^ Zq and mpopiii = 200Mq. The pre- 
dicted [O/H] at z = is in perfect agreement with the data. 
The scatter of the [Fe/H] distribution increases with z and 
becomes very large at z ~ 11 as a result of the large fluctu- 
ations in the merging history (see Fig. 1) and enrichment of 
halos among different realizations at z > 10. This result is 
consistent with the fact that Q = 1 at 2 = 11 (see Sec. 3.6) 
and therefore no unpolluted regions are expected below that 
redshift. The scattering induced by the stochastic nature of 
the merger tree process is similar, although different in na- 
ture, with respect to the abundance fluctuations produced 
by inhomogeneous mixing, a process not describe by our 
perfect mixing approximation. This similarity can be appre- 
ciated by comparing our results with those found by Mori 
& Umemura (2006) using a hybrid N-body/hydrodynamic 
code. These authors conclude that the metallicity is highly 
inhomogeneous during the early phases of galaxy formation 




5 10 15 20 

redshift 



Figure 4. GM elemental abundance as a function of redshift for 
the fiducial feedback model (see Table 1) and Zcr = W~'^Zq, 
^Poplll = 200Mq . Solid lines show the evolution of oxygen and 
iron abundance, averaged over 200 realizations of the merger tree; 
the shaded area delimits the ztlcr dispersion region. The [O/H] 
contribution by Pop III (dotted line) and Pop II stars (dashed 
line) are also shown separately. The point is the measured [O/H] 
in high- velocity clouds (Ganguly et al. 2005). 

[t < 0.1 Gyr for a IO^^Mq galaxy) and becomes more uni- 
form at lower redshifts. 

We can see that Pop III stars dominate the GM enrich- 
ment for 9 < z < 20 while Pop II stars are the dominant 
enrichment channel for < 2 < 9. The GM enrichment 
history is the result of 4 phases in which different physical 
mechanism dominate; these phases are marked in Figure 4. 
The first one (phase a), in the redshift range 17 2 ^ 20, is 
characterized by a very rapid increase of the elemental abun- 
dances. Given the small binding energy of the high-redshift, 
low-mass halos, both SN^y^, and the less energetic SNII ex- 
plosions can overcome the gravitational pull of parent halos, 
contributing to the GM enrichment. As the typical halo mass 
scale grows, as a result of hierarchical formation process, 
metals produced by SNII are retained inside galaxies and 
only SN^T- regulate the GM enrichment (phase b). However 
the abundance grow is limited by the fact that an increasing 
fraction of the GM is polluted above the Zcr, hence quench- 
ing the formation of Pop III stars which eventually comes 
to an end at the beginning of phase (c) when the whole GM 
is enriched to Z > Zcr- Thanks to the higher star formation 
efficiency of larger halos. Pop II stars, during this phase, can 
eject their heavy elements. Their ability to do so, though, 
decreases with time because of the counteracting effect of 
the gravitational potential of larger galaxies that are form- 
ing below 2 = 5 (phase d). Note that the GM bears the 
nucleosynthetic signatures of (predominantly) Pop HI stars 
at 2 > 10; their abundance pattern becomes more and more 
elusive at later times when Pop II stars start to dominate 
the enrichment. This physical picture remains qualitatively 
unaffected by changing the value of Zcr although the relative 
duration and amplitude of the different enrichment phases 
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Figure 5. Lower panel: Ratio of metals ejected by halos M^j as 
a function of their mass and redshift witli respect to the total 
amount of metals predicted in the GM at z = 0, M^°'. Curves 
represent Mej/M^?* = 5 X 10~(*'^'-^' isocontours; also shown are 
^min and Mmax- The Pop III stars termination epoch (2 = 11) 
is shown by the horizontal line. The two rectangles identify the 
position of the maxima (see text). Top panel: Cumulative frac- 
tional contribution, F^j , to M^?' integrated over redshift by halos 
with different mass (shaded area); points with associated itlcr er- 
ror bars represent stellar-to-total mass ratios in the corresponding 
mass bin (multiplied by 10). 

may vary. This enrichment history is particular important 
to interpret the observed MDF features as we will discuss in 
Sec. 6 as the various phases can be traced by the number of 
stars that are formed during each of them. 

5.3 Mass ejection from progenitor halos 

In Fig. 5 we show, for a single representative realization, 
the ratio of metals (or, equivalently, gas^^) ejected by halos, 
Mej, as a function of their mass and redshift with respect 
to the total amount of metals predicted in the GM at z = 0, 
M*°*. In the Figure, the curve Mmin {Mmax) denotes the 
minimum (maximum) halo mass in which star formation can 
develop, i.e. Mmin = M4,{z). Coherently with ACDM mod- 
els, at high redshifts (15 < 2 < 20) halos have typical masses 
close to Mmin. As the redshift decreases, more massive ha- 
los are produced (predominantly) via merging events which 
represent the dominant formation channel for z <Q, making 
Mmax ^ Mmin - As a general rule, the GM metal enrichment 
is dominated for z > 2 by low-mass galaxies; it is only at 
later epochs that the contribution of larger galaxies becomes 
important. In fact the distribution seen in Fig. 5 is bimodal; 
the first, more extended peak is found for 3 < 2 < 7 cor- 
responding to halos in the mass range M = 1O*~^M0; the 

tt These quantities are equivalent because of the assumption that 
the wind has the same metallicity as the ISM of the parent galaxy 
(see Sec. 3.5). 



second peak corresponds to larger galaxies M — 10^^ Mq 
at 2 < 2. The bi-modality can be understood as the re- 
sult of the larger aej of low-mass halos (due to their lower 
gravitational potential) and the larger gas content of more 
massive ones. The first effect dominates at high redshift; 
below 2 = 2 the ejected mass from larger halos becomes 
substantial despite their relatively low values of aej. Halos 
in the intermediate mass range are born from a relatively 
small number of merging events of gas-poor progenitors and 
therefore httle gas is left for further ejection. The largest 
halos instead are characterized by higher baryonic fractions 
as they incorporate a larger number of progenitors with dif- 
ferent histories. For example, halos with mass lO^-^Mo at 
2 = 5 have contributed alone to 0.5% of the metals present 
in the GM today. An equivalent contribution is provided by 
the most massive MW progenitors in the second peak. In the 
first peak a large plateau is presents around the maximum 
indicating that low-mass galaxies M < 6 x l(f Mq play a 
fundamental role for the enrichment of the GM. 

This is even more evident from an inspection of the 
top panel of the same Figure where we show the cumu- 
lative contribution to A/*°' integrated over redshift (Fej). 
The curve grows rapidly for 10* M© < M < 6 x 10® 
and M > 10^^ Mq, with a fiatter behavior in between, re- 
fiecting the afore-mentioned bimodal distribution. About 
60% of the ejected mass at 2 = comes from halos with 
M ^ 6 X 10® M0, compared to 20% contribution from the 
most massive galaxies. We conclude that GM enrichment 
is dominated by low-mass halos M ^ 6 x 10® M©. This re- 
sult is insensitive to the value of Zcr, although the relative 
contribution of Pop III and Pop II stars may vary. 

As in our code we also store the information about the 
stellar mass corresponding to each DM halo we can rephrase 
the previous results in term of such quantity (points in the 
top panel of Fig. 5). Proto-galaxies that mostly contribute 
to GM pollution have a typical stellar-to-total mass ratio 
A/. /A/ < 0.03, or Af* < 2 X 1O**A^0. 

Note however that the scatter of the Af*/Af ratio can 
be relatively large as shown by the error bars in the plot; the 
dispersion is caused by the different star formation histories. 
The scatter is maximum around M = 10® Af©; this popula- 
tion is the most numerous one through most of the MW 
history and they may be produced by widely different com- 
binations of formation processes and virialization epochs. 

From Fig. 4 we have concluded that Pop 111 dominate 
the enrichment down to 2 = 11; from Fig. 5 we can also 
set an upper limit to the fraction of today's GM metals 
provided by Pop III stars. At 2 = 11 the largest halo from 
which metals can escape is M < 2 x 10^ Mq, corresponding 
to F^^°^''''' < 0.04 i.e. Pop 111 stars contribute negligibly 
to the heavy elements currently detectable in the Galactic 
environment. 



6 DATA COMPARISON 

We are now ready to compare the results of our model with 
the best currently available data on halo metal-poor stars in 
terms of their MDF. Such distribution gives the number of 
MW relic stars as a function of their iron abundance, [Fe/H], 
a quantity commonly used as a metallicity tracer. In Fig. 6 
we compare the joint HK/HES MDF by Beers & Christlieb 
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Figure 6. Le/t panel: Comparison with the MDF for the galactic halo stars ob(served by Beers & Christlieb (2006) with the inclusion 
of the two hyper-metal poor stars (Christlieb et al. 2002; Frebel et al. 2005) (points) and the MDF obtained for the fiducial feedback 
model with Zcr = ZQ,mpopiii = 200Mq (histogram). The histogram is the average value of the MDF over 200 realizations of the 
merger tree and re-normalized to the number of observed stars with [Fc/H]< —2. The shaded area represents ztlu Poissonian errors. 
Right. The same results plotted in terms of the cumulative number of stars below a given [Fo/H]. 



(2006) with the simulated MDF (plotted both in differen- 
tial and cumulative form) obtained using our fiducial model 
with Zcr = W~'^Zq and mpopiii = 2OOM0. For complete- 
ness, we have added to the above sample the two known 
hyper-metal poor stars (Christlieb et al. 2002; Frebel et al. 
2005) using the new determination by Christlieb, Bessell & 
Eriksson (2006) for the iron-abundance of the star HE0107- 
5240. Due to contamination by disk stars, we have cut the 
data sample for [Fe/H]> —2. To compare with data, we have 
normalized the simulated MDF to the total number of ob- 
served stars with [Fe/H]< —2. Note that comparing these 
quantities we are implicitly making the assumption that all 
simulated metal-poor stars reside in the halo: this hypoth- 
esis is supported by the most recent N-body simulations 
(Scannapieco et al. 2006). 

The model reproduces the observed MDF quite well, 
particularly for [Fe/H]> —3.2. A marginally significant 
deviation from the data (always within 1-a) is seen at 
lower [Fe/H] i.e. a range populated by stars the formed in 
accretion-enriched halos. The two hyper-metal poor stars 
are not reproduced by our model. By definition in fact, the 
formation of Pop II stars is allowed only if the gas in the 
star-forming region has a metallicity Z ^ Zcr = 1Q~'^Zq 
adopted in this case. 

6.1 Changing the critical metallicity 

In order to bracket the current uncertainty on the determi- 
nation of Zcr, we have explored the sensitivity of our results 
to values of Zcr different from the reference Zcr ~ W~'^Zq 
(Fig. 7). Independently of the adopted Zcr value, the ob- 
served global properties of the MW are reproduced with the 
same accuracy as in the case Zcr ~ 1Q~'^Zq. We stress again 
that the choice Zcr = is equivalent to state that stars (in- 



cluding metal-free ones) form at all times according to a 
Larson IMF (see Section 3). 

The first striking result is that the MDF cutoff is shifted 
at gradually lower values of [Fe/H] when Zcr decreases (for 
reasons discussed in Sec. 5.1); the high-[Fe/H] limit of the 
distribution appears unchanged, showing, in all cases, a very 
good agreement with the data for [Fe/H]> —3.2. However, 
the tendency of the data to drop faster than than the model 
at lower [Fe/H], already noticed for Zcr = 10"^^©, persists. 
In particular, when Zcr = Zq the MDF shows a steep 
increase in —5 <[Fe/H]< —4.6 followed by a fiat plateau 
extending up to [Fe/H]= —4.1, a markedly different shape 
with respect to the previously analyzed case Zcr = 10~*Zq. 
This evolution bears the imprint of the CM chemical evolu- 
tion and to understand it in detail it is useful to recall the 
four enrichment phases identified in Fig. 4. During phases 
in which the CM metallicity grows slowly (as in phase b), 
stars with approximately equal [Fe/H] are formed; instead, 
when the CM metallicity changes rapidly (phases a and c), 
stars born at slightly different cosmological times from ac- 
creted gas may have quite different [Fe/H], giving rise to 
the MDF plateau. In general, though, the number of stars 
formed during phase (a) is too small to be appreciated from 
the normalized MDF. As Zcr is decreased from 10~*^q to 
zero, phase (b) contracts in time, being almost suppressed 
for Zcr = 0. The correspondingly lower CM metallicity, ex- 
plains the shift of the low- [Fe/H] MDF cutoff to lower val- 
ues. The metallicity increase during phase (c) brings [Fe/H] 
to values around -3, where Pop II always dominate the en- 
richment, independently of the value of Zcr- This physical 
picture holds qualitatively also for the lower Zcr case dis- 
cussed below. This result can hardly be overlooked, as it 
shows that the MDF embeds the features imprinted by me- 



© RAS, MNRAS 000, 1-17 



Cosmic stellar relics in the Galactic halo 13 



100 



10 




-4 -3 -6 -5 -4 -3 

[Fe/H] 

Figure 7. Same as left panel of Fig. 6 but for values of Zcr = IQ-^, IQ-^, Q Zq 



chanical feedback and reflects the heavy element enrichment 
history of the MW environment. 

As already mentioned, no stars have been detected so 
far by the various surveys in the range —5.3 < [Fe/H] < —4, 
a feature that we have dubbed as the "metallicity desert". 
Not only the model with Zcr ~ W~^Zq does not show such 
feature, but it also cannot account for the presence of the 
two HMP stars. The latter objects can instead be repro- 
duced if Zcr = W~^Zq. For the reasons explained above, 
the MDF plateau can extend now through the broad range 
—5.6 <[Fe/H]< —4.2. Even in this favorable case, though, 
the metallicity desert cannot be reproduced. If the desert is 
real or is a mirage produced by selection effects, it is im- 
possible to tell at this stage. What is clear though is that 
hierarchical models like the present one have serious diffi- 
culties to prevent the formation of stars in the desert [Fc/H] 
range. We will return to this in the discussion. 

Finally, the case Zcr = is somewhat peculiar. Al- 
though it predicts stars for [Fe/H]< —5.4, as expected from 
the previous arguments, its MDF shape is different with re- 
spect to cases with Zcr 7^ 0, showing a monotonic increase 
in the range —6 <[Fe/H]< —2. When Zcr ~ 0, in fact, no 
Pop III/II IMF transition occurs and, as a consequence, the 
GM enrichment proceeds more gradually. Even in this case, 
though, a considerable number of stars, is predicted in the 
metallicity desert. 



-2 


-6 


10-5, 10-' 


\OZq 


Zcr/ Zq 


Fo 


10'^Zq 


6.0 X IQ-'' 


10-6% 


5.8 X 10-^ 





7.5 X 10-3 



Table 2. Values of the _Fo parameter (see text) as a function of 
the assumed Zcr- 



Zcr = 10"" ^0, Zcr = 0. Models Zcr > 10"" ^0 give a value 
of _Fo compatible with the observational limit whereas for 
Zcr = we find Fo = 7.5 x 10"^ > FS'"' . In this case, in 
fact, the first burst of Pop III stars produces a huge amount 
of metal-free, long-lived stars which dramatically pushes _Fo 
beyond the allowed range. Thus, this parameter allows to 
robustly rule out models in which the critical metallicity 
is zero, and/or put an important constraint on primordial 
IMF. In the second option, in order not to conflict with the 
experimental limit, we must set a lower limit on the mass 
of the flrst stars equal to mpopiii > mi (2 ~ 20) ~ 0.9Mq, 
I.e. the stars have a mass larger than the turn-off mass at 
z ~ 20, the epoch at which the first star formation events 
occur in our model. 



6.2 The parameter Fo 

The useful parameter, Fo, has been originally introduced 
by Oey (2003) and later improved by Tumlinson (2006). Fo 
is defined as the number of metal-free stars divided by the 
total number of stars with [Fe/H]< —2.5. This parameter 
allows to quantify the implications of the persisting non- 
detection of metal-free stars in the Galactic halo. The ob- 
servational limit on the value of Fo is given by the inverse 
of the total number of stars with [Fe/H]< —2.5 observed 
into the galactic halo (where the numerator of the ratio is 
fixed equal to 1 to give a non- vanishing limit). Using the 
preliminary data by Beers and Christlieb (2006) we find 
that Fo°'"' = 1/1152 ~ 8.7 x 10"*. In Table 2 we present 
the values of Fo obtained using the models Zcr ~ 10"** Zq, 



6.3 Second generation stars 

An additional information that can be extracted from our 
study is the fraction of "second generation" (2G) stars, de- 
fined as the stars which have been enriched only through the 
nucleosynthetic products of Pop III stars. From an observa- 
tional point of view it is crucial to identify 2G stars, because, 
by studying their elemental composition we could trace back 
the properties of the first stars that polluted them. 

In Fig. 8 we identify the MDF regions populated by 2G 
stars and their fractional number for three values of Zcr - For 
Zcr > W~^Zq 2G stars populate both the lowest and the 
highest [Fe/H] tails of the MDF, i.e. they formed both in the 
accretion and self-enrichment modes. If Zcr ~ 0, on the con- 
trary, they can be produced only through self- enrichment. 
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Figure 8. Impact of second generation stars for tlie fiducial feedbaclc model Zcr = 10~*,10~®,0 Zq. In each panel, the highest 
histogram show the total MDF; the smaller histograms represent the MDFs for the 2'''*-generation stars. The percentages show the 
fraction of 2"'^-generation stars with respect to the total in the ranges in which they appear. 



This is because in this case the first (Type II) SN are not en- 
ergetic enough to expel their metals which remain confined 
within the parent halos. 



More specifically, the dominant formation channel of 2G 
stars is accretion when Zcr = W~^Zq and self-enrichment 
when Zcr = W~*Zq, refiecting the different Pop III/II tran- 
sition redshift. In the case Zcr = W~^Zq the transition oc- 
curs rapidly as a result of the lower amount of metals neces- 
sary to reach the critical threshold. Hence, 2G stars can be 
formed through self- enrichment only in the few halos host- 
ing Pop III stars; however, in unpolluted halos accreting gas 
from the GM, 2G stars can continue to form until the time 
at which feedback from Pop II stars raises the GM metallic- 
ity substantially. For the same reasons, the trend is reversed 
when Zcr = W~'^Zq. 



Along with their location on the [Fe/H] axis, the per- 
centages shown in Fig. 8 quantify the relative number of 2G 
stars with respect to the total number of stars in the iron- 
abundance ranges in which 2G stars are found. Two points 
are worth noting. First, the number of 2G stars formed by 
self-enrichment is negligible; even in the case Zcr ~ 0, their 
relative abundance is only 1.6% of the total population of 
stars with —2.5 <[Fe/H]< —2.3; this means that only 10 
of the 612 halo stars observed by Beers & Christlieb (2006) 
in the above [Fe/H] range could be truly 2G stars. Second, 
it is only in the metallicity bins close to the critical metal- 
licity value that the fraction of 2G stars increases substan- 
tially, thanks to the accretion mode of formation operating 
at those metallicities. However, given the very low number of 
extremely metal-poor stars available, their absolute number 
is likely to be close to zero. In practice, given the Beers & 
Christlieb (2006) sample, out of the 1150 stars observed in 
the range —4 <[Fe/H]< —2.5, from the results of our model, 
quantified in Table 3, we conclude that, independently on 
the Zcr value, the probability to detect 2G stars is virtually 
zero. 



Zcr/ Zq 


Number of 2G stars 


10-4 


1.3 


lo-f* 


0.3 





6 X 10-2 



Table 3. Number of 2G stars expected in the Beers & Christlieb 
(2006) sample as a function of the assumed Zcr- 

7 SUMMARY AND DISCUSSION 

Using the newly developed Monte Carlo code, GAMETE, 
we have followed the gradual build-up of the stellar popu- 
lations and metal enrichment of the Milky Way along its 
past hierarchical evolution within the concordance ACDM 
model. Adopting simple but physically motivated prescrip- 
tions for star formation and chemical/mechanical feedback 
we have first calibrated our model to the observed global 
properties of the MW. Once the model free parameters (the 
star formation, e*, and the wind, e„, efficiencies) have been 
fixed through such calibration, the results have been com- 
pared against the most recent experimental data (Beers & 
Christlieb 2006) on the halo stars Metallicity Distribution 
Function. We have also explored different values of the crit- 
ical metallicity, Zcr, which governs the transition from the 
Pop III to Pop II star formation mode (e.g. Schneider et al. 
2002). The main results are the following: 

(i) The global properties of the MW are well reproduced 
both using a closed-box model (e* = 0.5, — 0, i.e. no 
winds) and a model including mechanical feedback (e, = 0.7, 
En, = 0.2); however, winds are required to dilute metals in 
the Galactic Medium and hence reproduce the wide range 
of [Fe/H] values characterizing halo stars (see Fig. 3). 

(ii) Stars with [Fe/H]< —2.5 form in halos accreting GM 
gas enriched by earlier supernova explosions. 

(iii) The fiducial (feedback included) model for which 
Zcr ~ W~^Zq, mpopiii — 200Mq provides a very good fit 
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to the MDF, including the cutoff observed at [Fe/H]= —4. 
However: the existence of the two HMP stars with [Fe/H]< 
—5 cannot be explained within this model. 

(iv) Models with Z^r < 10~^Zq can account for the two 
HMP stars at the price of ovorpopulating the "metallicity 
desert", i.e. the range —5.3 < [Fe/H] < —4 in which no stars 
have been detected so far. 

(v) The current non-detection of metal-free stars, giving 
a value of the parameter Fq = 1/1152, robustly constrains 
either Zcr to be larger than zero or the masses of the first 
stars to be mpopiii > O.9M0. 

(vi) The statistical impact of truly second generation 
stars, i. e. stars forming out of gas polluted only by motal-froo 
stars, is negligible. For example, in the Beers & Christlieb 
(2006) sample only (at best) 1-2 stars could retain such 
metal-free nucleosynthetic imprint. 

(vii) Independently of Zcr, 60% of metals in the Galac- 
tic Medium arc ejected through winds by halos with masses 
M < 6 X IO^Mq, thus showing that low-mass halos are the 
dominant population contributing to cosmic metal enrich- 
ment. 

The previous results provide a self-consistent and coher- 
ent physical scenario for the formation and evolution of the 
Galaxy in a cosmological context. The model also explains 
in detail how the transition to Pop H stars occurred and 
clarifies by what mechanisms (mostly self-enrichment and 
gas accretion) the various regions of the MDF have been 
populated. 

In spite of these successes not all is clear and additional 
work will be necessary to understand some puzzling issues. 
No doubt that the most outstanding one is the unexplained 
existence of the "metallicity desert" . The problem in its sim- 
plest form can be stated as follows: in order to produce low- 
mass stars with [Fe/H] < —5, the value of Zcr must be 
decreased to similar levels (or below); as a consequence, an 
uncomfortably large number of stars is predicted by such 
a model in the desert. Assuming that the extremely low 
[Fe/H] of the two stars (which are very CNO-rich and have 
a total Z ^ Zcr) is not a result of some peculiar nucle- 
osynthetic history, ways to quench the formation of stars in 
the desert region must be devised. If the Pop HI star for- 
mation episodes indeed spread their metals into the GM to 
reach a metallicity levels of [Fe/H] w —5, they could also 
heat the gas to about Th ^ 10^ K. In such preheating sce- 
nario (Madau, Ferrara & Rcos 2001) gas accretion on viri- 
alizing halos with Tyir < Th would be prevented, whereas 
star formation and self-enrichment could continue in already 
collapsed structures, which would increase the GM metal- 
licity through their winds. It is only when dark matter ha- 
los of sufficiently large mass become numerous that accre- 
tion onto them could start again, now producing stars with 
[Fe/H]> —4. Translated in redshift, and guided by Fig. 4, 
the metallicity desert could translate in a gap of Az « 5 at 
2 > 10 in which accretion is halted by the preheating, an- 
other aspect of the mechanical feedback, which once again 
is found to play a crucial role in the problem. 

The Holy Graal of the various halo star surveys has un- 
doubtedly been the detection of a truly rnctal-frcc star. As 
these objects have so far escaped detection, growing atten- 
tion has been posed in the recent years on probing the nucle- 
osynthetic patterns of extremely metal-poor stars, with the 



hope that these objects arc born from gas polluted purely by 
(a single ?) metal- free SN. Unfortunately, our study shows 
that these second generation stars are no less elusive than 
their predecessors, having number frequencies that are so 
low to be zero even in the largest available sample used in the 
present analysis. These findings open a different perspective 
on the recent data by Cayrel et al. 2004 (later interpreted 
by Chiappini et al. 2005). Using a sample of 30 stars in the 
range —4.1 < [Fe/H] < 1.7 they measured the abundances 
of 17 elements from C to Zn, finding that these stars have 
an extremely low abundance ratio scattering. This is ob- 
viously contrary to what expected if these stars are direct 
descendants of the first gcncration(s) of stars formed after 
the Big Bang. Our results indicate that the stars used for 
the analysis axe almost certainly not second generation stars; 
instead, they represent a more mature population in which 
the abundance scatter has been smoothed out by several 
rc-proccssing cycles of the gas. 

A number of assumptions and approximations have also 
been done, mainly dictated by our persisting ignorance of 
some physical processes. These also need to be improved in 
future studies. The most critical one is probably the perfect 
mixing approximation which may acclcratc the extinction of 
metal free stars with respect to the epoch z = 11 predicted 
by our results. This problem is partially alleviated by the 
spread induced by the stochastic nature of merger histories 
which appears to be similar to that found by sophisticated 
numerical simulations of mixing in individual galaxies (Mori 
& Umemura, 2006). An additional key hypothesis concerns 
the IMF of Pop III stars which has been taken to be either 
a (5-function in the SN-^^ mass range MOAf© < mpopiii < 
26OM0 or a standard Larson IMF when Zcr = 0. Given 
the strong [Fe/H] dependence of mpopiii wc have assumed 
as reference value mpopiii = 200Mo but we have also ex- 
plored the implications of adopting the two extreme values 
of 14OM0 and 2601^© in the Appendix. A Gaussian Pop III 
IMF centered in 2OOM0, as suggested by some theoreti- 
cal works (see Nakamura & Umemura 2001; Scannapieco, 
Schneider & Ferrara 2003) might help explaining the con- 
comitant existence of the two HMP stars and the metallicity 
desert. 

Additional hypothesis have been made. Star formation 
has been take to occur in burst only, rather than in a contin- 
uous manner. This assumption is good for low-mass halos, 
but might become increasingly less accurate in larger galax- 
ies. However, the fact that the global properties of the MW 
are well reproduced, seems to indicate that this is not a 
major concern. Also, instantaneous recycling has been used. 
This is definitely a fair approximation for very short-lived 
Pop 111 stars, but might be a relatively poor description of 
the chemical evolution driven by more evolved generations. 
This might induce an overestimate of the metallicity, and 
consequently, an earlier Pop III/II transition than the actual 
one; when averaged over the present Hubble time, though, 
this problem should not affect the derived MDF properties. 

We finally like to comment on the differences between 
our results and those obtained by Tumlinson (2006). The 
treatment of chemical and mechanical feedback is probably 
more physical in our scheme, allowing metals/matter trans- 
fer between halos and the GM. Within our model, we find 
that the Pop III/II transition is a relatively sudden event, 
rather than a smooth one extending down to ^: < 6 as Tum- 
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linson found. This might be the result of the perfect mix- 
ing approximation that we have used, whose broad validity 
is supported by the afore-mentioned analysis of the metal 
filling factor evolution. Another difference between the two 
models resides in the calibration method: wo select the best 
models in terms of the two free parameters e* and as those 
reproducing simultaneously all the observed global proper- 
ties of the MW, rather than fixing them separately. As an 
indicator of the gas metallicity, we have used observations 
of gas in the Local Group (as for example the High Velocity 
Clouds), rather than the general IGM which is very likely 
not representative of the Galactic environment we are inves- 
tigating. 

In spite of these differences, there are no major tensions 
between the results of the two studies. It is worth noting 
that Tumlinson also identifies Zcr ~ W~*Zq as the most 
probable transition threshold, although his analysis could 
not discriminate between largely different Zcr values, given 
the poorer sample statistics; also, the two HMP stars are 
not taken into account in his analysis. 
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APPENDIX A: VARYING THE POP III MASS 

The models in Sec. 6 assume mpopiii = 200Mo and Zcr > 
0; we now explore the effects of setting mpopiii equal to 
140, 26OM0, respectively. These two values arc the extremes 
of the SN-y-y mass range and correspond to the minimum and 
maximum iron yield of Y^e = 2.8 x 10"" and Yfe = 0.45. 
We recall instead that Y = 0.45, independently of mpopiii- 

Following the same arguments using in Sec. 5.1 wo have 
checked that models without feedback largely fail to repro- 
duce the overall MDF shape also for mpopiii = 140, 26OM0; 
then need for feedback is independent of the adopted Pop III 
stellar mass. We have derived the MDF for galactic halo 
stars for these two cases; this is shown in Fig. 9 for the 
usual values of the two free parameters e, = 0.7, tw = 0.2. 

For mpopiii = I4OM0, the shape of the MDFs is com- 
pletely independent of Zcr and closely resembles the behav- 
ior found for the case Zcr = 0. This result from the extremely 
low Ypc of 14OAf0 stars which makes their contribution neg- 
ligible to that of Pop II stars which are then controlling the 
MDF shape. 

If mpopiii ~ 26OM0 instead, the MDF cutoff are quite 
sensitive to the value of Zcr, shifting to towards lower [Fe/H] 
values with decreasing Zcr ■ Such trend is very similar to the 
one already found for models with mpopiii = 2OOM0. How- 
ever in this case, given the higher Ype, the cutoff is located 
at higher [Fo/H]. For this reason the observed data arc bet- 
ter reproduced by model with Zcr < 10^^ ^0. Note however 
that even these models cannot account for the presence of 
the two HMP stars. 

In conclusion the MDF is sensitive to the value of 
mpopiii- However, as long as mpopiii is within the SN^-y 
mass range, it is always possible to find satisfactory fits to 
the data by suitable changing Zcr within its uncertainty in- 
terval Zcr = W'^'^^Zq. 

Most importantly both the value of Fo (Sec. 6.2) and 
the number of 2G stars (Sec. 6.3) are practically independent 
of the value of mpopiii- In particular, using the two best 
models mpopiii = 14OM0, Zcr = W~^Zq and mpopiii = 
26OM0, Zcr = W~^Zq we obtain Fo = 6.88 x 10"^ and 
6.3 X 10~®, respectively. The expected number of 2G stars in 
the Beers & Christlieb (2006) sample, to be compared with 
Tab. 3, is equal to 3.45 x 10"^ and 2.2 x 10"^ respectively. 
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Figure Al. Same as left panel of Fig. 6 but for different values of mpopiii = 140, 26OM0 and Zcr = IQ-*, 10"^, W'^Zq. 
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